arXiv:astro-ph/0503438vl 21 Mar 2005 


Astronomy &: Astrophysics manuscript no. 2629 February 2, 2008 

(DOI: will be inserted by hand later) 


Variation of the line profile moments for stars pulsating in 
distorted oblique non-radial modes 

O. Kochukhov 

Department of Astronomy and Space Physics, Uppsala University, SE-751 20, Uppsala, Sweden 
e-mail; oleg@astro.uu.se 

Received / Accepted 

Abstract. We derive expressions and develop a numerical technique for the analysis of the line proHle moment 
variations for stars pulsating in oblique non-radial modes. This represents an extension of the widely used spec¬ 
troscopic moment mode identification method to the oblique distorted pulsations observed in rapidly oscillating 
Ap stars. We demonstrate that a non-axisymmetric superposition of the pulsation and rotation velocity fields 
results in a qualitatively new behaviour of some of the line profile characteristics. It is found that for the majority 
of roAp stars the second moment varies with the pulsation frequency rather than with its first harmonic even 
for axisymmetric modes. We also identify pulsation observables which do not change during pulsation cycle but 
are modulated by the stellar rotation and can contribnte to the variability of the stellar spectra averaged over 
many pnlsation cycles. As an illustration of the new version of the moment technique, we compute rotational 
modulation of the pulsational changes of the line proHle moments for the obliqne axisymmetric dipolar pnlsation 
modes with different parameters. It is also shown that a distortion of the oblique dipolar modes predicted by 
the recent theoretical studies of the stellar magneto-acoustic oscillations can be readily diagnosed through the 
moment analysis. In particular, the shape of the pulsation phase modulation for the radial velocity and the second 
moment is very sensitive to non-axisymmetric pnlsation components, whereas the rotational modulation of the 
second moment amplitude is best suited to reveal axisymmetric magnetically induced distortion of pnlsations. 
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1. Introduction 

Analysis of the stellar pulsations is a valuable astrophysi- 
cal technique offering a unique possibility to constrain the 
internal structure and evolution of stars. A detailed in¬ 
formation about the surface geometry of the stellar pulsa¬ 
tions is a necessary prerequisite for a successful asteroseis- 
mic investigation into the physical properties of the stellar 
layers. Pulsational perturbations in the surface layers of 
stars are described by the spherical harmonic eigenfunc¬ 
tions. The angular degree £ and the azimuthal number m 
is sufficient to fully characterize the surface geometry of 
non-radial pulsations in normal slowly rotating stars. 

A range of the mode identification techniques was de¬ 
veloped to extract information about the pulsation geom¬ 
etry from the photometric 1 Watson 1 198811 and, most re¬ 
cently, spectroscopic (Aerts & Eyers 20001 Telting 1200811 
observations of the stellar pulsations. The moment tech¬ 
nique, introduced by Balona (I1986II and further developed 
by Aerts et al. and Briquet & Aerts II2()()8II . is one of 

the most widely applied and efficient spectroscopic mode 
identification methods. It attempts to diagnose the struc¬ 
ture of pulsational perturbations by interpretation of the 


time variation of the first three velocity moments that are 
linked to the line centroid (radial velocity), line width and 
skewness (line asymmetry). 

In contrast to normal non-radially pulsating stars, a 
considerably more complex and challenging problem is 
posed by the pulsations in rapidly oscillating Ap (roAp) 
stars (Kurtz & Martinez I2()n()ll . Since their discovery 
(Kurtz I1982II roAp stars have drawn great attention as 
the only example of the stellar oscillations in the presence 
of a global strong magnetic field. The stellar magnetism is 
closely linked to the excitation of the high-overtone short- 
period p-modes in roAp stars and, therefore, asteroseismic 
analysis may contribute to solving the long-standing prob¬ 
lem of the origin and structure of global stellar magnetic 
fields. 

The pulsational amplitude and phase of the oscil¬ 
lating magnetic stars were found to vary synchronously 
with the rotational modulation of the line of sight mag¬ 
netic field. A straightforward phenomenological explana¬ 
tion of this behaviour was given by the oblique pulsator 
model lKurtz n~982ll . which postulated that pulsations are 
aligned with the oblique magnetic field and modulation 
of the pulsational characteristics is a result of the pe- 
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riodic changes in the aspect angle at which pulsations 
are observed. In the standard oblique pulsator framework 
the modes were represented by the axisymmetric dipolar 
(£ = 1, m = 0) eigenfunctions. However, detailed photo¬ 
metric monitoring of some of roAp stars (e.g. Kurtz et al. 
nwi and theoretical studies (Shibahashi & Takata llbbdl 
Bigot & Dziembowski I2()(?^ Saio & Gautschv 120(1411 have 
suggested that pulsation modes in roAp stars are signifi¬ 
cantly distorted by rotation and magnetic field and can¬ 
not be described by a single spherical harmonic function. 
Recent time-resolved high-resolution spectroscopic stud¬ 
ies of roAp stars (e.g. Kanaan & Hatzes 119981 Kochukhov 
& Bvabchikova 12111113 also revealed a startling diversity 
of the pulsational behaviour of individual spectral lines, 
which was interpreted to be an interplay between a short 
vertical length-scale of pulsation waves and the presence 
of significant vertical gradients of chemical abundance 
in the line forming atmospheric regions (Kochukhov & 
H.va.bchikova ridllT^ Ryabchikova et a,1. 12 (1(1211 . 

These discoveries emphasized a necessity to develop 
quantitative tools to analyse in detail the surface geom¬ 
etry of pulsations in roAp stars in order to distinguish 
the horizontal and vertical pulsation effects. The measure¬ 
ments and modelling of the line profile moments appear 
especially well-suited for Ap stars which are typically very 
slow rotators. However, the latest version of the standard 
moment method implemented by Briquet & Aerts Il2()()dll 
cannot be directly applied to this kind of stars because 
this technique does not account for the mode obliquity 
and neglects magnetic distortion of pulsation modes. 

The main aim of the present investigation is to develop 
a numerical technique suitable for modelling the moment 
variations of stars pulsating in distorted oblique non-radial 
modes. Moreover, we characterize the main features of the 
pulsational variability of moments in the context of mod¬ 
ern theories of the magneto-acoustic stellar oscillations. 
The rest of the paper is organized as follows. Sect. [3 dis¬ 
cusses mathematical formalism and numerical implemen¬ 
tation of the extended moment technique. Calculations of 
the moment variations are presented in Sect. |3 for a vari¬ 
ety of oblique pulsation geometries and Sect. 0] concludes 
the paper with the summary and discussion. 

2. Influence of pulsation and rotation on the line 
profile moments 

In the present study we use representation of the pulsa¬ 
tion velocity suggested by Kochukhov Il20()4a|l . This pa¬ 
per demonstrated that an arbitrary pulsation velocity field 
corresponding to a mono-periodic linear stellar oscillation 
with the angular frequency w can be written as 

V{e, <j), t) = V^{6, (j>) cos {ojt) + V^{e, (j>) sin (wt), (1) 

where 6 and (f) are the spherical coordinates on the stel¬ 
lar surface, whereas V^{9, (j)) and V^{6, (p) are the vectors 
describing the surface distribution of the pulsation ampli¬ 
tudes for the horizontal and vertical velocity fluctuations: 

V\9, P) = {VP, V^} and V^{9, P) = {Vp, Vi, V^} .( 2 ) 


The complete expressions of the VH ^ velocity compo¬ 
nents in terms of the spherical harmonics were given by 
Kochukhov Il2004a|l . 

In the general case of a rotating star with oblique non- 
radial pulsations each oscillation mode is represented by 
a superposition of many spherical harmonic components 
(Bigot & Dziembowski I2()()2I Saio & Gautschv I2()()4II and 
there are no coordinate transformations which are able 
to reduce this description to a single spherical harmonic 
eigenfunction assumed in the classical moment technique 
of Balona llT?Mll and Aerts et al. (nnniii- As a consequence, 
the six pulsation velocity components V^g^ have to be de¬ 
fined in any coordinate system through the sums over the 
spherical harmonic functions Y^m(^, </>) and their deriva¬ 
tives. Nevertheless, there still exists a special reference 
frame, e.g. the one connected with the global magnetic 
field of the star, where the pulsational displacement ex¬ 
hibits certain symmetries and the respective mathemati¬ 
cal description of the velocity field is the least complicated 
(Saio & Gautschv l2()()4|) . We choose this coordinate sys¬ 
tem, inclined by the angle (3 relative to the stellar rotation 
axis, to define the local angular coordinates 9 and (p and to 
specify pulsation amplitude maps V^{9,(j)) and V^{9,(j)) 
in terms of the sums over the spherical harmonic compo¬ 
nents. 

For the purpose of numerical computation of the disk- 
averaged characteristics of the line of sight velocity the 
visible stellar hemisphere is divided into a number of sur¬ 
face zones with approximately equal areas (as in Piskunov 
& Kochukhov Position of each surface element in 

the stellar reference frame is characterized by the latitude 
p (—7r/2 < p < 7 r/ 2 ) and the longitude I (0 < I < 27r). 

The transformation of an arbitrary vector field from 
the spherical coordinate system not aligned with the stel¬ 
lar rotation axis to the Descartes reference frame of the 
observer was treated by Piskunov & Kochukhov in 

the context of modelling the stellar magnetic field. Similar 
expressions are applicable for the pulsation velocity field. 

The transformation of the vectors and ’P’* can be writ¬ 
ten as a series of matrix multiplications employing the 
operator T ( 7 , S) defined as 

/ cos 5— sin^ 0 \ /l 0 0 \ 

T( 7 ,i 5 )= sind coS(5 0x0 cosy —siny , (3) 

yO 0 1J yOsiny cos y J 

where the angle y describes a tilt of the polar (z) axis of a 
given coordinate system and the angle S corresponds to a 
rotation around this axis. Thus, vector components of the 
pulsation velocity amplitudes in the observer’s coordinate 
system whose z-axis is directed along the line of sight are 
evaluated as 

/ v,^’%p,i)\ ( v;’\9,p)\ 

Vi’^ip,l) =T{-i,0)T{P,p)T{9--,ct>) -Vi’^{9,p) .(4) 

Here ip = fit is the phase angle corresponding to uniform 
stellar rotation with the angular frequency fl and i is the 
angle between the stellar rotation axis and the line of sight. 
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The total time-dependent velocity with respect to the 
observer includes the contributions from pulsation and ro¬ 
tation: 

K {p,l) = {p, 1) cos (uit) + V/ {p, 1) sin (ujt) + (p, 0 ■ (5) 

The 0 component of the rotation velocity is given by 

V^°^{p, 1) = cos p sin {I + (p)ve sini, ( 6 ) 

where Ve sini is the projected rotational velocity. 

We use the usual assumption that the intrinsic line 
profile can be approximated by a Gaussian with variance 
. This approximation will fail only for a (relatively un¬ 
common) situation of a magnetically sensitive line in the 
spectrum of a strongly magnetic and very slowly rotating 
roAp star. We follow the moment dehnitions of Aerts et 
al. and express variability of moments in terms of 

changes with pulsation frequency w, as well as its first and 
second harmonics, 2lo and 3w. The expressions for the ra¬ 
dial velocity (V), the second moment {V"^) and the third 
moment (V^) thus become 

{V) = (t4(p, Z)) = Vioocos(wt) -f Void sin (wt), (7) 


with the angular frequency uj and rotation with the angu¬ 
lar frequency Q. The latter is included in the expressions 
indirectly through the time dependence of the coordinate 
transformation defined by Eq. 0 . Taking into account 
this fundamental background of the problem of investiga¬ 
tion into the nature of oblique non-radial pulsators it is 
natural to take an advantage of the stellar rotation and 
to observe the pulsation velocity field at different aspect 
angles - a unique possibility realized only in roAp stars. 
Then the mode identification for roAp stars and the anal¬ 
ysis of the distortion of pulsations by rotation and mag¬ 
netic field can be carried out by modelling the temporal 
evolution of the pulsation spectra of (V), (V^) and (V^) 
as a function of the stellar rotation phase, similar to the 
studies of the rapid photometric variability of roAp stars 
(Kurtz & Martinez I2()()( III . To describe the theoretical ro¬ 
tational modulation of each amplitude t4(v^) in Eqs. 0- 
0 we compute an equivalent representation in terms of 
the total time-dependent amplitude Ay^{ip) and the phase 
^’x(V5) 

Kx(t) = cos {nivt)+ b^{ip) sin {nu)t) , , 

= A^{(p) cos [nujt + $x(</5)], 


(V^) = cr^-I-(K/(p, 0) = + M^2 oo + 1^020 + 2 K 002 ) 

-I- 2 V 101 cos (wt)-I-2Voii sin (wt) 



-b 

\ (E200 

— ^020) COS ( 2 u;t) + V110 sin 

{2ujt) 

and 





(E 3 ) 

= 

(K3(p,Z))+3a2(K(p,0) = |(E 2 oi 

+ E021) 


-b 

I (E300 

+ AV 102 + E120 + 4 cr^E 100) 

0 

0 

T 


-b 

I (Eo 30 

+ 4 E 012 + E210 + 4 ct^Eoio) 

sin [iot) 


-b 

1 (^201 

— E021) cos {2u)t) -b 3 Eiii sin {2u}t) 


-b 

\ (E300 

— 3 Ei 2 o) cos ( 3 a;t) 



-b 

\ ( 3 E210 ~ Eo3o) sin ( 3 tut). 



(9) 


The quantities Vknm represent the surface integrals of 
various powers of the V^, and velocity components 
and are approximated with the weighted sums over the 
visible surface elements 

N 

VuuM = E 0]" [vHp, 1)T [v^ip, ny ■ (10) 

i=l 


Here the weight function Wi is evaluated under the as¬ 
sumption of the linear limb-darkening law 


IE, = 


3 (1 - M -I- upi) piSj 
7 r (3 — u) 


( 11 ) 


where u is the limb-darkening coefficient, whereas pi and 
Si are the cosine between the surface normal and the line 
of sight and the intrinsic area of the Ah surface zone re¬ 
spectively. In deriving expressions 00 we took into ac¬ 
count the symmetry of the rotation velocity field which 
reduces the Eqoj sums to zero for all odd values of j. 

We emphasize that the moment equations 00 in¬ 
clude variability on the two different timescales: pulsation 


where 

^ = yjal{ip) + bl{ip) and $x(7’) = •(^^) 

Examination of the moment formulae derived for the 
oblique pulsators reveals some essential differences with 
respect to the moment behaviour in normal pulsating 
stars. Apart from the expected rotational modulation of 
all terms in the moment amplitudes and phases, the most 
important feature is the relative amplitudes of the terms 
contributing to the pulsational variation of the second mo¬ 
ment. In stars where pulsations are aligned with the rota¬ 
tion axis and the pulsation velocity is of the same order 
as the rotation one, contribution of the main frequency 
term is usually small or comparable to the 2uj term and 
even reduces to zero for any rotation velocity if pulsa¬ 
tions are axisymmetric (i.e., in the language of terms of 
Eq. 0: Eioi,oii(/3 = 0, m = 0) = 0). This does not hap¬ 
pen for oblique pulsators. Eor all but the slowest rotat¬ 
ing roAp stars pulsation amplitudes are smaller than the 
projected rotational velocity. Hence, Hioi.oii(/3 7 ^ 0) » 
(E 200 — E 020 ) even for an axisymmetric pulsation and we 
expect to see a strong single wave variability dominating 
the pulsation behaviour of (E^) in most of the roAp stars. 

Another interesting feature of the extended moment 
equations is the rotational modulation of the “constant” 
terms in the expression for the second moment and the 
presence of such terms for the third moment. Although 
changes of these terms on the timescale of the stellar ro¬ 
tation are relatively small, this effect illustrates nicely a 
consequence of the non-asymmetric co-addition of the pul¬ 
sation and rotation velocity fields and may contribute to 
the variation of line prohles in the regular (averaged over 
many pulsation cycles) spectroscopic observations of roAp 
stars. 
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In comparison with the classical semi-analytical ex¬ 
pressions of the moment variation derived by Aerts et al. 
(nnnsi, the numerical scheme introduced here is consid¬ 
erably more general and straightforward to implement. 
It is also more flexible due to the possibility to include 
additional physical effects (such as stellar surface inhomo¬ 
geneities) in the weight function W. A similar numerical 
approach was also followed in the recent modification of 
the moment method presented by Briquet & Aerts 
who optimized the moment technique for multiperiodic ro¬ 
tating stars but did not address the problem of modelling 
the moment variation for stars with oblique pulsations. 

3. Illustration of the moment variations 

For all calculations of the line profile moments presented 
in this section we used a 1876-element surface grid, lin¬ 
ear limb-darkening coefficient u = 0.5 and adopted an 
intrinsic profile with cr = 5 kms“^. The maximum sur¬ 
face amplitude of the pulsation velocity was set to Vmax = 
max ^= 4 kms“^. These values corre¬ 
spond to a typical limb-darkening, thermal Doppler line 
width and pulsational amplitude of the late-A rapidly os¬ 
cillating peculiar stars. The reference values of the other 
relevant parameters were chosen to be close to those of 
the prototype roAp star HR 3831 (Kochukhov et al. 120041 
Kochukhov I2004bl) : Ue sinf =30kms“^, i = 60°, (3 = 90°. 
We also attempt to test an effect of a significant horizon¬ 
tal pulsation component by comparing calculations with 
the ratio of the horizontal to vertical pulsation amplitude 
AT = 0.5 (cf. Saio & Gautschv r2004ll with the case of purely 
vertical pulsations. 

Figures HHS present rotational modulation of the am¬ 
plitude and phase of each term in the variation of the first 
three moments. 

3.1. Axisymmetric dipolar modes 

The classical oblique pulsator model is represented by 
a pure axisymmetric dipolar (^ = 1, m = 0) mode in¬ 
clined with respect to the rotation axis iKurtz 1198211 . If 
the mode obliquity is close to 90°, the radial velocity vari¬ 
ation reaches maximum at the two phases during rotation 
cycle (Fig. corresponding to the moments of the best 
visibility of each of the two pulsation poles. The amplitude 
modulation is clearly non-sinusoidal, as it features sharp 
minima and broad maxima. 

In agreement with the qualitative discussion in Sect. [3 
in the case of a moderately rapid rotation {ve sini ^ Fmax) 
variation of (F^) is dominated by the main frequency term 
and occurs with the tt/2 radian shift in rotation phase rela¬ 
tive to the modulation of radial velocity. As demonstrated 
by Fig. n the variation with twice the pulsation frequency 
becomes dominant only when the maximum pulsation ve¬ 
locity exceeds UeSin*. Modulation of the respective term 
in the amplitude of the second moment is sinusoidal and 
occurs out of phase with respect to modulation of the main 


frequency term. Harmonic contributions to the pulsational 
changes of the third moment are always small. 

A significant horizontal pulsation displacement with 
the ratio of the horizontal to vertical pulsation amplitude 
K — 0.5 has no effect on the shape of the rotational mod¬ 
ulation of the amplitudes of moments. At the same time, 
horizontal motions have a profound effect on the maxi¬ 
mum amplitude achieved during rotation cycle: the (F^) 
amplitude decreases by roughly a factor of 2, whereas the 
maximum amplitude of the third moment increases by the 
same amount. 

For an intermediate (/3 = 60°) obliquity of the axisym¬ 
metric dipolar mode we preferentially sample one of the 
pulsation poles. Consequently, one of the maxima of the 
amplitude of (F) and (F^) decreases. The rotational mod¬ 
ulation of the amplitude of the second moment is weakly 
affected. 

The pulsational phase modulation of the axisymmetric 
oblique dipolar modes is illustrated in the first and second 
columns of Fig. |3 The appearance of the phase modula¬ 
tion curves does not depend on Ve sin i and the presence of 
the horizontal pulsation fluctuations. The phase of radial 
velocity exhibits two tt radian phase jumps characteris¬ 
tic of the oblique dipolar pulsation. The rapid change in 
phase occurs when pulsation equator is passing through 
the plane defined by the rotation axis and the line of sight. 
In contrast, the tt radian jump in the phase of the main 
frequency term of the second moment occurs when pulsa¬ 
tion poles are facing the observer. The phase of the first 
harmonic term is not modulated by the stellar rotation, 
whereas the respective term in the pulsational changes of 
the third moment shows four jumps during rotation cycle. 

An effect of the moderate mode obliquity on the phase 
modulation is illustrated in the second column of Fig. 13 A 
decrease in the visibility of one of the pulsation poles does 
not change the overall character of the phase modulation, 
but modifies the times of the tt radian jumps, so the latter 
are no longer equidistant in rotation phase. 

3.2. Non-axisymmetric dipolar modes 

In their recent theoretical study Bigot & Dziembowski 
ii^nia have introduced a modified oblique pulsator model 
of roAp stars, which represented a significant departure 
from the classical axisymmetric dipolar pulsation geom¬ 
etry. Bigot & Dziembowski suggested that an ac¬ 

curate non-perturbative treatment of the interaction be¬ 
tween stellar rotation, pulsations and magnetic field may 
lead to a representation of the pulsation eigenmodes with a 
complex superposition of the spherical harmonic functions 
which contain substantial non-axisymmetric components. 
Furthermore, strong effects due to the centrifugal force 
are expected to break locking of the pulsation axis with 
the dipolar magnetic field, so that pulsations are aligned 
with neither the magnetic or rotation axes. However, this 
entirely new geometrical picture of the roAp oscillations 
remains somewhat speculative, as Bigot & Dziembowski 
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Fig. 1. Rotational modulation of the pulsation amplitudes of the radial velocity (V), second moment (V^) and third 
moment (F^) computed for the oblique dipolar (i = l,m = 0) non-radial mode with Fmax = 4 kms“^ and i = 60°. 
The other relevant parameters - the projected rotation velocity Ve sini, the ratio of the horizontal to vertical pulsation 
amplitude K and the mode obliquity /3 - are indicated at the top of the uppermost panel in each column of plots. 
Different curves show amplitudes for the terms changing with pulsation frequency (solid line), its first (dashed line) 
and second (dotted line) harmonics. 


Il2()()2ll gave no general predictions for the amplitude of 
the non-axisymmetric contribution to the pulsation veloc¬ 
ity field and failed to consider distortion of the pulsation 
modes due to a strong magnetic field typical for the ma¬ 
jority of roAp stars. 

It turns out that analysis of the moment variation pro¬ 
vides a very sensitive diagnostic of a non-axisymmetric 
pulsation. In this section we investigate an effect of small 
non-axisymmetric (m = ±1) contributions to the dom¬ 
inant oblique dipolar mode. In the pulsation reference 
frame the amplitude of the non-axisymmetric components 
is set to 10% of the velocity due to the £ = 1, m = 0 
component. This distortion of the dipolar pulsation has 
a marginal effect on the modulation of the amplitudes of 
moments (Fig.|21 first column). Rotational changes of the 
radial velocity amplitude become smoother and the am¬ 
plitude does not drop to zero at phases 0.25 and 0.75. At 


the same time, even very small non-axisymmetric com¬ 
ponents have a dramatic effect on the phase modulation 
curves, as illustrated in the third column of Fig. 0 The tt 
radian jumps vanish and are replaced by the steep but con¬ 
tinuous phase changes. This phenomenon is a straightfor¬ 
ward consequence of the surface distribution of the pulsa¬ 
tion velocity phase. All axisymmetric modes have constant 
phase with the tt radian jumps at the node lines, whereas 
the non-axisymmetric pulsation components are charac¬ 
terized by a continuous change of the pulsation phase from 
one surface zone to the next. This is readily reflected in 
the disk averaged velocity moments and allows to iden¬ 
tify a very low amplitude non-axisymmetric component of 
oblique pulsations. 




















6 


O. Kochukhov: Variation of the line profile moments for oblique non-radial pulsators 


;=l.m=-lAl.K=0.0,^=9tf' /=l^.m=0,K=0.0,^=90” /=l^,m=0, K=0A^=60" /=1,3, m=0, K=0.5.;3=90f’ 



ao oj (u 0.6 i.a ao 0.2 04 0.6 0^ 1.0 oo a2 0^ as ixt a2 u os 0.8 ix* 

Rolatiao phase Rocatioa phase Rotatioa phase lUitatioa phase 


Fig. 2. Same as Fig. ^ for the amplitude modulation of the moments computed for distorted dipolar modes. First 
column: non-axisymmetric dipolar mode with Vm=i = Vm=-i = 0.1V5„=o, second to fourth columns: magnetically 
distorted axisymmetric mode with 14=3 = 0.514=i and different K and j3. 


3.3. Magnetically distorted dipolar modes 


In another important theoretical development Saio & 
Gautschy have studied the effect of a strong mag¬ 

netic field on the geometry of oblique non-radial modes. 
The stellar rotation and the resulting non-axisymmetric 
pulsational contributions were neglected, but calculations 
were performed for the magnetic field strength of up to 
9 kG. In the oblique pulsator model of Saio & Gautschy 
<1200411 pulsations remain axisymmetric and are aligned 
with the magnetic field axis. The effect of dipolar mag¬ 
netic field is to constrain oscillations to the surface areas 
close to magnetic poles. Thus, an initially dipolar pulsa¬ 
tion is distorted in such a way that it is described by a 
sum of the odd axisymmetric pulsation components. We 
computed moment variation for such a pulsation model 
assuming that the octupolar, £ = 3, pulsation has a 50% 
smaller amplitude compared to the £ = 1 component and 
contribution of other harmonics is negligible. 


The amplitude modulation of the line profile moments 
computed for magnetically distorted oblique pulsations is 
presented in Fig. |21 (columns 2-4). The most significant 
effect of the £ = 3 component appears in the modulation 
of the second moment amplitude. The respective curve be¬ 
comes more complex and two additional minima develop 
at the rotation phases 0.25 and 0.75. Variation of the sec¬ 
ond moment is also quite sensitive to the obliquity of the 
mode and to the presence of the horizontal pulsation mo¬ 
tions (predicted in the theory of Saio & Gautschy I2004II . 
Thus, magnetic distortion and related enhancement of the 
horizontal amplitude of pulsations in roAp stars can be 
most efficiently studied with the analysis of the rotational 
modulation of the variation of the second moment of spec¬ 
tral line profiles. 

Despite important influence on the moment ampli¬ 
tudes, a moderate contribution of the octupolar axisym¬ 
metric spherical harmonic considered here does not change 
modulation of the phase curves (Fig. 13 fourth column): 
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Fig. 3. Rotational modulation of the pulsation phase of the line profile moments. The format of the figure is similar 
to that of Figs. Hand 12 


it remains identical to the phase modulation obtained for 
the pure oblique dipolar mode. 

4. Summary and discussion 

Recent advances in theoretical studies of the stellar 
magneto-acoustic pulsations emphasize a necessity to ob¬ 
tain detailed empirical information about the structure of 
the pulsational disturbances in the atmospheres of roAp 
stars. Until recently, analyses of this interesting class of 
non-radial pulsators were limited to the time-resolved pho¬ 
tometric observations. Interpretation of the broad-band 
photometric variability turned out to be a challenging 
task because of the difficulty to discriminate between 
various effects contributing to the luminosity variation, 
the anomalous wavelength dependence of the photomet¬ 
ric pulsation amplitudes and phases fWatson ll988ll . non¬ 
linearity and non-adiabatic effects (Medupe 1200211 . As a 
result of these complications, little progress has been made 
in understanding the vertical and horizontal geometry of 
the roAp pulsations beyond the generic oblique pulsator 


model of Kurtz CnH2. In fact, the recent study by Saio & 
Gautschy (Enni has demonstrated that, due to the aver¬ 
aging over the visible stellar hemisphere, the photometric 
observables are very weakly sensitive to the details of the 
pulsation geometry and turn out to be practically useless 
in testing predictions of elaborate theoretical models. 

The newly emerged field of the observation and inter¬ 
pretation of the pulsational line profile variations in roAp 
stars (Kochukhov & Ryabchikova 12001 a^ offers a more 
promising research direction. Due to their simple relation 
to the pulsational displacement, the spectroscopic observ¬ 
ables provide a more direct, reliable and general approach 
to modelling the roAp pulsations compared to interpreta¬ 
tion of the photometric pulsational data, which requires 
some restrictive assumptions about the physics of pulsa¬ 
tions to be made. The most detailed and straightforward 
way to reveal the pulsation geometry of roAp stars is to 
construct a two-dimensional image of the pulsation veloc¬ 
ity field from the line profile variation using the pulsation 
Doppler imaging technique (Kochukhov l2(K)4al I2004bll . 
However, this sophisticated modelling requires spectro- 
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scopic observational data of a superb quality and is limited 
to moderately rotating roAp stars. An alternative possibil¬ 
ity is to apply the moment method (Briquet & Aerts l^nUHIl 
which is based on interpretation of the low-order moments 
of the absorption line profiles and can be applied to very 
slowly rotating stars. 

In this paper we presented an investigation of the diag¬ 
nostic potential of the moment method in its application 
to the oblique distorted stellar pulsations. We extended 
the moment technique to deal with the stellar pulsation 
velocity field with the symmetry axis other than the axis of 
stellar rotation and developed a numerical method for the 
computation of the rotational modulation of amplitudes 
and phases of different terms in the variation of the line 
profile moments. Calculations with our novel version of 
the moment technique reveal some important differences 
in the behaviour of the line profile moments in normal and 
oblique pulsators. In particular, for all but the slowest ro¬ 
tating roAp stars with slightly distorted dipolar modes we 
expect to observe a single wave variation of (V^) (the sec¬ 
ond moment) during pulsation cycle^. Furthermore, a non- 
axisymmetric superposition of the pulsation and rotation 
velocity distributions results in new terms in the variation 
of the higher line profile moments which are modulated 
by the stellar rotation but are constant on the pulsation 
timescale. 

Calculation of the line profile moments for various pul¬ 
sation geometries of roAp stars predicted by the theoret¬ 
ical studies demonstrates that analysis of the rotational 
modulation of variablity of the first two moments ((I^) and 
{V^)) is sufficient to diagnose various non-axisymmetric 
and axisymmetric distortions of the basic oblique dipolar 
pulsation geometry and allows to carry out an extremely 
sensitive verification of the theories attempting to describe 
an interaction between pulsation, magnetic field and stel¬ 
lar rotation. For instance, the non-axisymmetric modified 
oblique pulsator theory of Bigot & Dziembowski II2()()2II 
can be very efficiently tested for a large number of roAp 
stars by using observations of the rotational modulation of 
the pulsation phase of radial velocity and higher line pro¬ 
file moments. The presence of a smooth phase modulation 
instead of the tt radian jumps is a unique signature of non- 
axisymmetric modes. It reveals a weak non-axisymmetric 
contribution even in the background of a much stronger 
axisymmetric component. Similarly, the presence of the 
odd axisymmetric harmonic components arising due to 
the influence of the global dipolar magnetic field (Saio & 
Gautschv 120041) can be verified by the observation of the 
second moment variation. 

The promising diagnostic potential of the new moment 
technique can be realised through interpretation of the 
high resolution observations of the line profile variation in 
oblique pulsators. The bright, high-amplitude roAp stars 

^ The presence of such single wave second moment variation 
in the extremely slowly rotating roAp star 7 Equ (Kochukhov 
& Rvabchikova l 2 ()()lat is thus confirmed to be a signature of a 
non-axisymmetric non-dipolar pulsation. 


with short rotation periods (^10'^) are the best candidates 
for the application of the moment method because their 
surface pulsation velocity structure can be observed from 
different aspect angles during a few night observing run. 
Observational data for the prototype roAp star HR 3831 
have already been acquired by Kochukhov & Ryabchikova 
Il2nnibll . The time-series analysis of this object and in¬ 
terpretation of the pulsational spectral variation with the 
moment technique will be presented in forthcoming paper. 
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